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Abstract. 

A comparative study between two distinct perturbative series expansions for 
the pair-creation contact process is presented. In contrast to the ordinary contact 
process, whose supercritical series expansions provide accurate estimates for its critical 
behavior, the supercritical approach does not work properly when applied to the pair- 
creation process. To circumvent this problem a procedure is introduced in which one- 
site creation is added to the pair-creation. An alternative method is the generation of 
subcritical series expansions which works even for the case of the pure pair-creation 
process. Differently from the supercritical case, the subcritical series yields estimates 
that are compatible with numerical simulations. 
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1. Introduction 

The critical properties of nonequilibrium systems have attracted attention and have seen 
a great development in the last years [H [2|, [3] . This may be credited to their resemblance 
with critical phenomena occurring in systems in thermodynamic equilibrium. In 
equilibrium as well as in nonequilibrium phase transitions there is a singular dependence 
of the steady-state properties upon the control parameters, marking a transition between 
distinct regimes. Other common features include long-range correlations, well-defined 
order parameter and singularities characterized by critical exponents that are grouped 
in universality classes [4j. 

Due to absence of a general theory for the nonequilibrium regime, these properties 
are analyzed in particular models, trying to form a complete picture. In particular, 
a class of simple models used to perform this task is constituted by lattice models 
with absorbing states. The presence of an absorbing state is sufficient for breaking 
the detailed balance condition, making these models intrinsically irreversible. In the 
steady-state regime, these models display a phase transition from an absorbing state to 
an active state, as the control parameter is varied. 

A very robust universality class for these models is the directed percolation (DP) 
class. It describes the critical properties of any system with a phase transition between 
an absorbing and an active state characterized by a scalar order parameter, short-range 
interactions and no conservation laws [5]. 

Usually, the tool utilized in the study of these problems is the numerical simulation, 
but analytic approaches are also used as an important tool in this task leading, in some 
cases, to very precise results. As in other fields of physics, simulational and analytical 
approaches are complementary in the study of phase transitions in equilibrium and 
nonequilibrium systems [6]. One of these anal5d;ical approaches is the series expansion 
[7j that has provided the most precise estimates of the DP exponents in one-dimension 
[8]. Inspired by the equilibrium case, the series for nonequilibrium systems can be either 
supercritical or subcritical, in analogy with the high- and low-temperature expansions 
in equilibrium systems. 

A pioneer work that uses series expansions for lattice models with absorbing states 
was developed by Dickman and Jensen [7j, in their study of the critical properties of the 
contact process (CP). The contact process is an icon among the models that display a 
phase transition between an absorbing and active state. Initially the CP was proposed 
as a model for spreading of an epidemic disease [9] and has become the "Ising model" 
for the DP universality class. Using subcritical and supercritical series, Dickman and 
Jensen obtained the critical point and its associate exponents. Their results indicate 
that the supercritical case works better than the subcritical series. 

Recently, de Oliveira [lOj proposed an alternative approach for generating the 
subcritical expansion obtaining critical values comparable to those of the supercritical 
series [7]. The difference between the approaches is that in [10], the non-perturbated 
operator (associated with the annihilation of particles) should be diagonalized, whereas 
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this procedure is not required for the supercritical series. 

In this paper, a comparative study between subcritical and supercritical series 
expansion is performed for the pair-creation contact process (PCCP) \TT\. Although 
this model does not show any surprise concerning the critical properties, since it belongs 
to the DP universality class as the ordinary CP, it is interesting for the comparative 
analysis. We will show that the original approach for the supercritical series [7j does 
not work properly, presenting problems in its construction. Even an alternative version 
of the model, that would overcome these problems, does not lead to precise results for 
the critical point. On the other hand, the subcritical series gives results, without any 
maneuver, that are comparable to those obtained by numerical simulations [12] . 

This paper is organized as follows. In section 2 we present the model and the 
operator formalism. In section 3 we show the derivation of the supercritical series, 
discussing its problems. In section 4, we attempt to overcome the problem by introducing 
a modified PCCP and deriving the supercritical series. In section 5 we present the 
approach for the subcritical case and results coming from the Pade analysis for both 
pure and modified PCCP. Finally, in section 6 the conclusions and final discussions of 
this work may be found. 

2. Pair-creation contact process and operator formalism 

We consider an interacting particle system on a one-dimensional lattice with L sites. The 
system evolves in the time according to a Markovian process with local and irreversible 
rules. The configurations are described in terms of occupation variables, rji, with 
rji = 0,1 according to whether the site i is empty or occupied by a particle, respectively. 
The time evolution of the probability P{r], t) of a given configuration rj = (r/i, 772, rji) 
at the time t is given by the master equation, 

|p(r/, t) = j:^{w,{^^)P{r,\ t) - w^Piv, t)}, (1) 

where rj^ = 772, 1 — rji, rji). The transition rate Wi{ri) for the PCCP is given by 

Wiiv) = ^{^- Vi){Vi-2Vi-i + Vi+iVi+2) + Oir]i. (2) 

If an empty site has two occupied neighbor sites, then this site becomes occupied with 
a transition rate equals to n/2, where n is the number of the pairs of the first neighbors 
occupied. If the site is occupied, then it is becomes empty with a rate a, independently 
of its neighborhood. 

When the lattice is entirely empty, the system is trapped in an absorbing state. 
However, for small values of the rate a, an active state can be achieved in the stationary 
state. In this way, a continuous phase transition occurs in this model with the critical 
point localized, according to simulational results [12], in ac = 0.1341(6). 
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To develop this operator formalism, we represent the microscopic configurations of 
the lattice by the direct product of vectors 

\v) = <S>h)- (3) 

i 

The algebra is defined by the creation and annihilation operators for the site i 

A\Vi) = r]i\r]i - 1). (4) 
In this formalism, the state of the system at time t may be represent as 

\m) = j:piv,m- (5) 

{n} 

If we define the projection onto all possible states as 

(l=E(^l' (6) 
iv} 

the normalization of the state may be expressed as ( | ■?/') = 1. 

The master equation ([1]) can be shown to be equivalent to the following time 
evolution equation 

j^\m) = s\m), (7) 

where S is the evolution operator, given by 

S = Sp + aSa, (8) 

where 

Sp = l E{(4 - a4) (4-2^.-24-1 A-i + 4+iA+i4+2A,+2)} (9) 

is the pair-creation operator and 

Sa = E(^. - 4A) (10) 

i 

is the annihilation operator. 

3. Supercritical series expansion 

By rescaling the time, it is possible to rewrite the evolution operator as 

S = V + iJ,W, (11) 

where W = Sa, V = 2Sp and n = 2a. Since the operator W is associated to annihilation 
process and the creation of particle is present in the operator V , for small values of the 
parameter the creation process is favored, so that the composition above is convenient 
for a supercritical series. The action of each operator, in a general configuration (C), is 
explicitly shown by the expressions 

W^(C) = E(CO-r(C), (12) 

i=l 
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pi P2 

V{C) = + 2E(Cr) - {Pi + 2p,){C), (13) 

i=l i=l 

where r is the occupied sites number, (C-) is obtained replacing a particle at the site i 
by a hole, while pi and p2 are the number of the empty sites with one and two occupied 
pair of neighbors, respectively. Finally, (Cf) is the configuration obtained replacing the 
hole at the site z by a particle. 

Since we are interested in the steady-states properties, it is convenient to take the 
Laplace transform of the state vector, given by 

ms)) = / e-'^m))- (14) 

Jo 

Inserting the formal solution given by the equation (I7l), we find 

ms)) = {s~S)-'\m)- (15) 
The stationary state lip^oo)) = lim^^oo |^(^)) then be found noticing that 

|7/;(oo)) = lims|7/i(s)), (16) 

which is obtained integrating equation (fMl) by parts. Assuming that {ipis)) may be 
expanded in powers of and using equation f|T5l) . we have that 

ms)) = iV^o) + /ilV^i) + /i'|V^2) + --- = is-V- fiW)-'\m)- (17) 
Expanding the operator {s — V — fiW)~^ in powers of fi, 

{s-v-i2wy^ = (s-v)-^ [i + /i(s - vy^w + fi\s - vy^w'^ + ■■■], (i; 

and comparing each order in /i, we arrive at the expressions: 

|^o) = (^-l^)-^|^(0)), (19) 

and 

\^n) = {s-V)-'W\tPn-l), (20) 

for n ^ 0. 

The action of the operator [s — V)~^ on an arbitrary configuration (C) may be 
found noticing that 

(. - V)-\C) = S-' {(C) + {s- V)-'V{C)} , (21) 

and using the expression (fT3l) for the action of the operator V, we get 



[s-V)-\C) = sJ{C) + {s-V)-' 



Pl P2 

E(cn+2E(c;) 

i=i j=i 



(22) 



where Sg = l/{s + q) and q = pi + 2p2. 

The main obstacle to generate the standard supercritical series for this model is due 
the action of the operator (s — over a configuration (V) without pairs of particles, 
whose result is given by 

{s-V)-\V) = l{V). (23) 
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The successive action of {s — V)~^ over the configuration {V) or their off'springs will result 
in new configurations whose associated coefficient will be proportional to {n > 1) 
and thus, it diverges in the stationary state. To overcome this difficulty, we present in 
the next section an alternative version of PCCP model that avoid this divergence and, 
in principle, enables us to find asymptotic results for the critical properties of the PCCP 
model. 



4. Pair-creation contact process modified by a single creation transition 



The modification consists of adding to the transition rate ([2]) a term of creation by single 
particles similar to the ordinary contact process [1], so that the transition rate is given 
by 

u^iiv) = ^(1 - Vi){Vi-2Vi-i + Vi+iVi+2 + liVi-1 + Vi+i)} + (^Vi- (24) 

The operator V in equation (fTTj) now reads V = 2{Sp + •jSs), where the single creation 
operator Sg is given by 



Ss = l E{(4 - aaI){aUa_, + 4+1 A+i)}- 



(25) 



In the limit of 7 — 0, we recover the pure PCCP model. 

The operator W acts over (C) according to equation f|T2l) . whereas the action of the 
operator {s — V)~^ is given by 



{s-V)-\C) = s,,,{{C) + {s-V) 



-1 



pi 



P2 



i=l 



i=l 



+ 



+ Sp,^-f\is-V)-^ 



li=l 



i=l 



(26) 



where Sp_^ = (s + [pi + p2\ + 7[ri + ^2])"^ pi and p2 (are the number of empty sites 
with one and two pairs of occupied neighbors, respectively and ri and r2 are the number 
of empty sites with one and two occupied neighbors, respectively. The expression {C'-) 
corresponds to the configuration obtained by replacing a hole at site i by a particle. 
Differently from the pure PCCP, the modified PCCP does not have divergences when 
one takes the limit s — 0. Thus, for configurations (P), we have that 



1 



76 



(27) 



where h is an integer number and {V^) is a configuration originated from (V) by a simple 
creation process. The particular case of the vacuum (0) for which (s — V)^^(O) = l/s(0) 
does not diverge in the stationary state limit. 

In the supercritical expansion the operator (s — V)~^ acting over any configuration 
(except the vacuum) generating an infinite set of configurations, so it is impossible 
to evaluate \il^{s)) in a closed form. We can evaluate, however, the ultimate survival 
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probability Poo, which corresponds to the coefficient of the vacuum. We define the series 
coefficients as 



Using a proper computational algorithm, we can calculate the coefficients up to order 
23, with values of 7 shown in figure [H The limiting factor for this calculation is actually 
the memory required. We attempt to obtain the critical properties of the PCCP model 
by performing the numerical extrapolation 7 — 0, since this approach does not work in 
the case 7 = 0. 

To analyze the series, we use the d-log Pade approximants approach. These 
approximants are defined as ratios of two polynomials 



In our case the function /(a) represents the series for ^ InPoo(a), for a fixed value for 
the parameter 7. Therefore, we are able to obtain approximants satisfying the condition 
L + M < 22. One verifies that diagonal {L = M) and near-diagonal approximants 
usually exhibit better convergence properties, so that we will restrict our calculations 
to the set of approximants such that L = M + 6, with 6 = 0, ±1. 

In the neighborhood of the critical point the ultimate survival probability behaves 
as Poo ~ (ttc — a)^- The critical point ac is determined by the pole of the approximant 
Flm{c() and the exponent P by the residue associated with this pole. 

In figure [H the estimates for critical points ac are depicted for different values of 
the parameter 7. Using a linear extrapolation of these data, we obtain the estimation 
for the asymptotic limit 7 ^ as ac = 0.1398(7), in disagreement with the simulational 
result etc = 0.1341(6) [12]. We remark that decreasing the parameter 7, the dispersion 
of the approximants increases and the precision of the result becomes smaller. 

The exponent f3, as shown in figure [2], is in accordance with the DP universality 
class value for any value of the parameter 7, as expected. 

However, the uncertainties on the location of the critical point put in doubt the 
accuracy of the supercritical series around 7^0. In fact, what happens is that the 
coefficients of the configurations (V) in the steady-state behaves as 



diverging as 7 ^ 0. This behavior causes a ill-conditioned series for 7 0, becoming 
impossible a more accurate result for the critical point in this limit. 

5. Subcritical series expansion 



Poo = 1 - 5I«i(7)a'- 





/(«)• 



(29) 




(30) 



To develop subcritical series expansion, we rewrite the evolution operator as 

S = W + X{Sj, + -fSs), (31) 
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Figure 1. Critical line in the plane a versus 7 for the modified PCCP model by 
supercritical series (circles). The dashed line is a linear extrapolation toward 7 = 0. 
For comparison, we plotted the numerical simulation result (cross). 
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Figure 2. Exponent (3 for the modified PCCP model as a function of 7. 



where A = a~^. Since the operator W is associated to annihilation process and, for small 
values of the parameter A, this annihilation process is favored, this expansion describes 
indeed the subcritical regime. Differently from the supercritical case, the series here are 
generated directly, without the necessity of taking the Laplace transform. 

The operator W, associated to annihilation particles process, is expressed as 
W ^ J2i Wi, with Wi ^ Ai - A\Ai. Each term Wj has the following set of right 
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and left eigenvectors 

|0) = |o), (o|^(o| + (.|, (32) 
with eigenvalue Aq = and 

|l)^-|o) + |.), = (33) 

with eigenvalue Ai = —1. 

To find the steady state vector \ip), that satisfies the steady condition Slip) = 0, we 
assume that 

oo 

e=i 

where \i/jo) is the steady solution of the non-interacting term W satisfying the stationary 
condition 

W\iJo) = 0. (35) 

The vectors {ipe) can be generated recursively from the initial state {ipo)- Following 
Dickman [13], we get the recursion relation 

|7/>,) = -i?(S'p + 75s)|^£-i). (36) 

The operator R is the inverse of W in the subspace of vectors with eigenvalues and is 
given by 

^= E \<f^n)^m, (37) 
n(^0) 

where |0„) and are right and left eigenvectors of W, respectively, with nonzero 
eigenvalue A„. Since the creation of particles is catalytic, if we start from steady state 
of the noninteracting term, that corresponds to the vacuum state, we will obtain a trivial 
steady vector. To overcome this problem, it is necessary to introduce a modification 
on the rules of the model. The necessity of changing the initial state in systems with 
absorbing states in order to get nontrivial steady states have been considered previously 
by Jensen and Dickman [71 [T3] and more recently by de Oliveira [10] . 

The modification we have performed [TU] consists in introducing a spontaneous 
creation of particles. For the case 7 = (pure PCCP) the creation occurs in two 
adjacent sites, chosen to be z = and i = 1. This modification leads to the following 
expression to the operator W 

W = Y.W, + q{Uo + U^-Wo- W,), (38) 

i 

where q is supposed to be a small parameter and Ui = A- — AiA\. The steady state l^/'o) 
of W is not the vacuum state anymore. Now, it is given by 

|V'o) = |.0.)+2g|.10.) + g2|.ll.), (39) 

where all sites before and after the symbol "." are empty. 

Two remarks are in order. First, only the last term in will give nonzero 
contributions to the expansion so that I'lpi), ^ > 1, will be of the order q^. Second, 
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Although the change in W will cause a change in R, only the terms of zero order in the 
expansion in q, given by the right-hand side of equation (|37|) . will be necessary since the 
corrections in R will contribute to terms of order larger than g^. For instance, the two 
first vectors, and \4'2), for the pure PCCP are given by 

l^i) = q'{2\.l.) + + 1.101.) + ^1-111.)}, (40) 

and 

1^2) = g'{^|.i.) + ^|.ii.) + ^l-ioi.) + ^|.iii.)+ 

+ ^1.1001.) + ^1.1101.) + ^1.1011.) + ^1.1111.)}, (41) 

where the translational invariance of the system is assumed. 

For 7 7^ 0, it suffices to consider the spontaneous creation in just one site. This 
modifications leads to the following expression 

W = Y.W, + q{Uo-Wo), (42) 

i 

instead of equation (155]) . 

From the series expansion of the state vector it is possible to determine several 
quantities. In this paper, we will be concerned only with the series expansion for the 
total number of particles N, given by 

N={.0.\J2n,\^). (43) 

i 

One can show that the coefficient of in the expansion for is simply the coefficient 
of |.l.) in \ipe). This allows us to get a longer series for the number of particles. For 
both modified and pure PCCP, we have obtained a series with 38 terms. 
The critical behavior of obeys the following relation [10] 

AT- (a, -a)-'' 1 1 (44) 

where and rj are the exponents related to the time correlation length and to the 
growth of the number of particles, respectively. Through an analysis of the d-log Pade 
approximants, like that performed for supercritical series expansion, we are able to 
calculate the critical point and its associated exponent. Using the series for the number 
of particles A^ the values obtained were ac = 0.1310(14) and i/||(l + rj) = 2.68(15). 
These values do not agree with the best estimates, ac = 0.1341(6) [12] and z/||(l + rj) = 
2.2777305 [8J. Another estimate is obtained by performing a non-homogeneous d-log 
Pade approximant ^5\. Using this approach we obtained the value ac = 0.1337(17), 
closer to the value obtained by numerical simulations [T2] . 

More reliable estimates are obtained when we determine the Pade approximants for 
the series {ac — a){d/da) InN = 9 [3116], called biased Pade approximants. Considering 
a given Pade approximant [L/M] and a trial value of etc, we develop the series above 
obtaining 9{a^. We can build curves for different Pade approximants by repeating this 
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Figure 3. Biased estimates oi 9 = i/||(l + 77) as a function of dc derived from Pade 
approximants to the series {dc — a){d/da) \nN = 9 evaluated at dc for the PCCP. 
The approximants shown are [17/18], [18/17], [18/18], [18/19], and [19/18]. 



procedure for several trials ac and we expect that they intercept at the critical point 

From the figure [3] we see a very narrow intersection of the Pade approximants, 
revealing the utility of this approach. However, as pointed out by Guttmann [17], it is 
difficult to estimate uncertainties in series calculations. Thus, in order to give a more 
realistic estimate of the quantities measured here and their associated uncertainties, we 
have estimated them by taking into account the first and last crossings among various 
Pade approximants. The estimate of ac for the PCCP, ac = 0.13396(1), in excellent 
agreement with the corresponding value obtained from recent numerical simulations 
|12j . Finally, the exponent obtained is z/||(l + r]) = 2.274(3), in agreement with the best 
estimation of z/||(l + r]) for the DP universality class [8]. 

To stress the difference between subcritical and supercritical approaches, we have 
also developed subcritical series expansion for the modified PCCP model for some 
values of 7, as shown in figure HI In the limit 7 — 0, the subcritical approach gives 
ac = 0.1340(1), which is in good agreement with the estimate ac = 0.13396(1) for 7 = 
and ac = 0.1341(6) obtained from numerical simulations [12] . On the other hand, the 
extrapolation for 7 — > in the supercritical case does not lead to the correct value, as 
can be seen in figure [TJ 
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Figure 4. Comparative results for tlie modified PCCP model using supercritical 
(circles), subcritical (triangles) series and numerical simulations (cross). 



6. Conclusions 

In this work, we have considered the pair-creation contact process by carrying out a 
comparative study between subcritical and supercritical expansions. Differently from 
the contact process, the supercritical expansion in its canonical formulation [7] is 
intrinsically ill-conditioned for the PCCP. In fact, divergences in coefficients of the 
supercritical series occur for any model with creation by cluster of particles. As an 
attempt to circumvent the divergences, we introduced and analyzed a modified version 
of the PCCP. However this procedure does not lead to the correct result when one 
recovers the original model. It would be interesting to search an alternative procedure 
to develop supercritical series for such models. We have shown here that subcritical 
series expansion provides results comparable with recent numerical simulations [12j. 

A natural extension of this work is the use of subcritical series to study models with 
creation by clusters (pairs and triplets) and diffusion of particles [HI |12] . Results from 
numerical simulations [12], [18] suggest that these models display a tricritical point. Thus 
this approach, associated with an analysis by means of partial differential approximants 
in two- variables, could be used to determine the existence of a tricritical point. In fact, 
such an approach has already been shown to be useful for the location of a multicritical 
point in a generalized contact process p^. 
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